An interacting resonant level coupled to a Luttinger liquid: Universality of 

thermodynamical properties 
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We investigate a model of a single resonant level coupled to the edge of a quantum wire in the 
Luttinger liquid phase or to the middle of a chiral Luttinger liquid via both tunneling and contact 
interaction. Utilizing the Yuval- Anderson approach, we map our model onto a classical ID Coulomb 
gas in which all the details of both the interactions in the lead and the level-lead interaction enter 
only through the corresponding Fermi-edge singularity exponent, which we explicitly evaluate using 
the Bethe ansatz solution for a particular model of the lead. Thus the population, dynamical 
capacitance and level entropy are universal (i.e., equal for different models), as we demonstrate to 
hold quantitatively using density matrix renormalization group calculations. Since the Coulomb 
gas description is of the single channel Kondo type, we infer that the universality we found implies 
that Luttinger liquid physics has no qualitative effect on these properties, contrary to perturbative 
results. 

PACS numbers: 71.10.Pm, 72.10.Fk, 73.21.Hb, 73.40.Gk 
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Introduction. — Understanding the properties of 
strongly correlated systems has been one of the main sub- 
jects of investigation in condensed matter physics in re- 
cent years. An important class of such problems is that 
of quantum impurities, i.e., systems with a finite num- 
ber of degrees of freedom coupled to reservoirs of non- 
interacting particles, the best known examples of which 
are the Kondo and Anderson models [lj and the spin- 
boson model 0. Another important type, for which 
non-Fermi liquid physics is well established, is that of 
one-dimensional electronic systems. When no symmetry 
is spontaneously broken, the low energy physics of those 
systems is described by the Luttinger liquid (LL) theory, 
where the quasiparticles are bosonic modes of density 
(or, in a dual description, phase) fluctuations It is 
then natural to try to bring these two themes together, 
by studying quantum impurities coupled to LLs: from 
the quantum impurity perspective, the reservoir now has 
a non-trivial physics of its own; from the LL point of 
view, this gives a way to probe the intricate physics of 
the electrons which are coupled to the impurity (and not 
the relatively simple behavior of the weakly-interacting 
bosons). Besides these fundamental motivations, quan- 
tum impurities (e.g., quantum dots) and LLs (e.g., quan- 
tum wires), are the basic ingredients of nano scale cir- 
cuits, so that understanding them has a profound im- 
portance for applications. Although such models have 
been studied theoretically for some time, most effort has 
been concentrated at understanding transport properties 
[H 13 j with a recent revival of interest in the application 
of exact methods to non-equilibrium situations [B|, @, 0] ; 
other phenomena have usually received only scant atten- 
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Model. — The simplest possible system to study these 
effects is that of a single level coupled to the edge of a 
LL (which can be realized, e.g., by lithographically defin- 
ing a small quantum dot at the end of a quantum wire, 



or by coupling an impurity atom to the end of a car- 
bon nanotube), or, equivalently, to the middle of a chiral 
LL (e.g., an anti-dot near the edge of a fractional quan- 
tum hall bar) [l3| . We include contact interaction be- 
tween the level and the lead. Here transport properties 
are not relevant; however, many other interesting ques- 
tions can be investigated. In this letter we concentrate 
on thermodynamical properties: The level population, its 
dynamical capacitance (which can be probed experimen- 
tally by capacitively coupling the system to a quantum 
point contact), and the level contribution to the entropy 
and specific heat. We demonstrate, both analytically and 
numerically, that these properties show universality, and 
depend on the different interactions only through a sin- 
gle parameter, the Fermi edge singularity exponent of the 
system. Using an analogy to the single channel Kondo 
problem, we show that this universality means that these 
physical quantities have an essentially Fermi liquid like 
behavior, which is affected only quantitatively, but not 
qualitatively, by the LL physics. 

The system is described by the following Hamiltonian: 



H = H icad {^(x),i>{x)} + e dU - [711^(0) + H.C.] 
+ [V f (0)V(0)- ^(0)^(0)], (1) 



where d, ip{x) are Fermi operators of the level and 
the lead, respectively, Sq is the bare level energy, 711 is 
the level-lead tunneling matrix element, and An is the 
strength of the level-lead interaction. 

Yuval- Anderson approach. — Using the Yuval- 
Anderson approach [l4|, in either the canonical [l5| or 
the path- integral [HI formulations, we expand the parti- 
tion function to all orders in the hopping clement in the 
Hamiltonian, and evaluate the resulting many-particle 
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Green functions. We then obtain the expression: 
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where £o is a short time (ultraviolet) cutoff, To is the 
(renormalized) level width (whose expression is given be- 
low), and (3 is the inverse temperature of the original 
problem. This expression thus has the form of a classical 
grand canonical partition function of a one dimensional 
gas of particles residing on a circle of circumference /3, 
with fugacity ^/Tq^o/it. Each particle is assigned a pos- 
itive (negative) charge if it corresponds to hopping of 
an electron from the lead to the level (vice- versa) . The 
charges must thus be alternating, with an overall charge 
neutrality. Thus, a configuration is completely specified 
by the sign of the first charge (denoted by a in the above 
expression) and by the positions of the particles. The 
action of this classical system consists of two terms: 

S({r<},*)= (-l) i+j a FES m{ , } 



+ e \/3^+a £ (-l)Sl. (3) 

[ l<i<2JV J 

The first term is an interaction between the particles, 
with the form of a Coulomb interaction between charged 
rods, and a coefficient (charge squared) ofeS) the Fermi 
edge singularity exponent of our problem (by which we 
refer to the scaling dimension of for vanishing hop- 

ping term). If the wire has a finite length L, the argu- 



at zero temper- 



ment of the logarithm is sin[7rVs{T] _ Ti)/L] 
ature (v s is the velocity of the LL's bosonic excitations), 
whereas at finite temperature it is replaced by an elliptic 
function [3j . The second term in the action of the classi- 
cal system corresponds to the energetic cost of filling the 
level, and resembles the effect of an electric field applied 
on the charges. 

We now discuss in more details the parameters appear- 
ing in the Coulomb gas model. In the case where the 
lead is noninteracting, we have the usual resonant level 
model, where it is known that «fes = (1 - H) and 
To = 7r |tu| 2 fo cos(5), where 5 — tan~ 1 (7ri'oAu/2) is the 
phase shift of the electrons in the lead caused by their 
interaction with the level, and i>q is the local density of 
states in the lead 
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When there are nonvanish- 
ing interactions both in the lead and between the level 
and the lead, the situation is more complicated. From 
bosonization Q we obtain (since there is no backscat- 
tering in this problem) ckfes = (1 — g\\/'Kv s ) 2 /g, and 



TABLE I: Parameters appearing in the Coulomb gas model, 
Eqs. © and (J3J). See the text for further details. 

Non-interacting lead Bosonization Suggested generalization 



o?fes (1 - §<5) 

To 7r|7u| ^ocos(<5) 



i(l-^l) 3 lh- g l5 eS ) 



7r|7ii| i>o 



t|tu| focos(5 c ff) 



To = 7r|7n| vo- If we take the limit of noninteracting 
lead, we see that in the bosonization treatment expres- 
sions that should contain the phase shift 5 are replaced by 
their leading order dependence on An. Thus, the values of 
apES an d r in any realistic model are renormalized by ir- 
relevant operators not appearing in the Luttinger model. 
On the other hand, «fes may be found, in general, using 
boundary conformal field theory arguments which 
relate it to the finite size spectrum of the lead with dif- 
ferent potentials at its edges. For a lead represented by 
a tight binding chain with hopping matrix element t and 
nearest- neighbor interaction U, finding these energies is 
equivalent, via the Jordan- Wigner transformation [3j, to 
finding the energies of an XXZ spin chain with boundary 
magnetic fields, which is exactly solvable b y th e Bethe 
ansatz [20| . Using the methods of Refs. |2jj, |22j, an ana- 



lytic expression for q;fes can be found in this case, which 
reads: 
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where g — 7r/[2cos _1 (— U/2t)], while t\\ and U\\ are the 
lattice analogues of level-lead hopping and interaction 
parameters of Eq. ([1]), to which they are related by the 
lattice spacing a: 711 = tny/a, and An = Una. It seems 

Un 



natural to identify <5 c ff — tan 



as an ef- 



y (2t)2-iP_ 

fective phase shift, which reduces to the usual phase 
shift when the lead is noninteracting. We may thus 
expect that for general models we can write «fes = 
I (l — ir^cff) i f° r some effective phase shift S e s (re- 
stricted to the range [— 7r/2, 7r/2]), so that Tq will be 
given by ir \ jn\ v o cos(<5 c ff ). In the following we will con- 
firm these results quantitatively by our numerical data. 
This discussion is summarized in Table Q] 

Universality. — We thus see that the partition func- 
tion depends on the original model only through three 
parameters: To, so and q;fes- The latter, in particular, 
contains the main effects of the interactions, both in the 
lead and between the level and the lead. This implies a 
universality in this system: many of its properties depend 
only on these three parameters, so that they will be the 
same for very different systems, with different strengths 
and signs of interactions, provided these three parame- 
ters are indeed the same. The properties which exhibit 
universality are the thermodynamical ones, e.g.: the level 




FIG. 1: (Color online) A color map of the differential ca- 
pacitance (divided by the noninteracting value) obtained by 
DMRG at half filling, as a function of the interactions in the 
lead and between the level and the lead, with contours of con- 
stant Q!fes superimposed. It can be seen that the effect of all 
the interactions comes only through this parameter, confirm- 
ing universality. See the text for further details. 

population and its correlation functions (or, equivalently, 
the static and dynamic level capacitance), and the level 
contributions to the entropy and the specific heat. 

We note in passing that the mapping into the Coulomb 
gas can be easily extended to include the case of an 
ohmic environment coupled to the level. The only ef- 
fect of this on the analysis is modifying the parameter 
c*fes by adding to it the impedance of the environment 
divided by the quantum resistance h/e 2 [§]. Hence, all 
our results apply to this case too. The universality is 
thus seen to have an even broader scope of applicability. 

Although the mapping to a Coulomb-gas applies rigor- 
ously only to the low frequency (long time) behavior, we 
have found it to hold quantitatively, at least when irrele- 
vant perturbations of the LL are not too strong. To this 
end we have performed density matrix renormalization 
group (DMRG) [23[ calculations on tight binding real- 
izations of the system discussed, with nearest neighbor 
hopping and interaction along the chain (of strengths t 
and U, respectively, as described above), keeping its pop- 
ulation at half filling. In Fig. [1] we show the differential 
capacitance dn/deo at £o = in a color map as a function 
of both the level-lead interaction and that in the lead. 
In all cases we have kept To = lCP 4 i and L = 50v s /t, 
modifying in and L accordingly, so as to keep all the pa- 
rameters of the Coulomb gas constant except «fes- On 
the color map we superimposed a contour plot of ofeSj 
taken from Eq. (j4|) . Indeed it is seen that the contours of 
constant ttFES are also contours of constant differential 
capacitance, confirming the important role of the former 
in determining the behavior of the system. Deviations 
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FIG. 2: (Color online) Level population as a function of its en- 
ergy: different symbols denote four models used in the DMRG 
calculations, while the different curves (which are a guide to 
the eye) correspond to different ofes values (the smaller qfes 
the wider the curve is and vice versa). In the last model 
the strengths of the nearest neighbor and next nearest neigh- 
bor interactions (in units of the hopping strength t along the 
chain) are: {1.5, 0.5}, {1.0, 0.5}, {0.5, 0.5}, {-0.25, -0.25}, 
{-0.5, -0.5}, {-0.75, -0.5}, in order of decreasing q:fes. In 
the third model U was taken as ±0.5£, with opposite sign to 
the corresponding fourth model case. See the text for further 
details. 

are seen only for quite strong interactions, where irrele- 
vant terms in the Hamiltonian are initially quite strong 
(and are not renormalized to zero because of the finite 
system size), and modify the results quantitatively. To 
appreciate this one should remember that for \U\ > 2t 
the system is no longer a LL [but becomes charge den- 
sity wave (phase separated) for positive (negative) £/]; 
whereas for \U\\\ > It the potential of ±[/n/2 felt at the 
last site of the lead when the level is full (empty) is strong 
enough to form a bound state. Both of these effects are 
not included in our treatment. 

A more detailed comparison is made in Fig. [5] Here 
we show the full dependence of the level population on 
its energy. The population curves corresponding to dif- 
ferent aFES values are presented, and on each such curve 
there are symbols of four types, denoting the numerical 
results on four different models: (a) A non- interacting 
lead with nonzero level-lead interaction; (b) A lead with 
nearest-neighbor interactions but zero level-lead term; (c) 
A system with both nonzero U\\ and nearest-neighbor in- 
teraction in the lead [which serves as a test to Eq. (01 and 
the subsequent discussion]; (d) A lead with next-nearest 
neighbor interactions (but vanishing U\\), used to show 
that our results apply even to non-integrable models (In 
this system g was determined numerically). The param- 
eters of the four models were chosen to give the same 
q:fes value [i.e., in each case we have chosen arbitrarily 
the interactions in the lead in models (c) and (d), and 
determined by the above condition all the other interac- 
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tions. The other parameters are the same as in Fig. [TJ 
except for the lead length, which is twice as large here]. 
Again we can see that the population is universal, de- 
termined by o:fes alone, and not by the parameters of a 
specific model. 

Lessons from the Kondo effect. — We now discuss an- 
other implication of the Coulomb gas mapping. The 
Coulomb gas we have obtained is similar to the original 
one, derived by Yuval and Anderson in their treatment 
of the anisotropic single-channel Kondo model In 
particular, the level population (minus one half) in our 
system is equivalent to the magnetization of the Kondo 
impurity, the level energy £o is analogous to a local mag- 
netic field, Tq plays the role of J±, and cufes is deter- 
mined by J z . We can thus immediately import all the 
known results from the Kondo problem [l[ to the case 
of a LL lead coupled to a level. The system considered 
can be in one of two phases - a strong coupling (anti- 
ferromagnetic Kondo like) delocalized phase, and a weak 
coupling (ferromagnetic Kondo like) localized phase. At 
very small values of Tq the transition is at cufes = 2, 
whereas for larger Tq it occurs for larger values of Ofes- 
In the localized phase, the low energy physics is that of 
an effectively disconnected level, so that its population 
is discontinuous as a function of Eq, and there is a non- 
vanishing residual entropy at zero temperature. Similar 
results regarding this phase, as well as the phase transi- 
tion line, were already discussed in Ref. [§], albeit using 
different techniques, and tested numerically by the au- 
thors 12j. On the other hand, in the delocalized phase, 
the impurity is well hybridized with the conduction band, 
so the level population is analytic in eq. One can write an 
explicit expression for the latter using the Bethe ansatz 
solution of the Kondo problem [l[ . In particular, for small 
values of Eq, one has: 



/ \ 1 £ Q 



(5) 



with Tk (the effective level width) corresponding to the 
"Kondo" temperature of the problem, which, for small 
r , is given by: 



Tk ~ (r £o) 



l/(2-a PES ) 



(6) 



and thus reduces to Tq for vanishing interactions 
(<2fes=1)- Hence, in this phase the population does not 
show any power law dependence on Eq. The only power 
law appearing is in the formula for Tk- However, the 
power depends on q;feSj an d is nontrivial (i.e., different 
from unity) even for a Fermi liquid lead if level-lead in- 
teractions are not negligible. In addition, in this phase 
the correlation function of the level population as a func- 
tion of time (denoted by r) will decay as (tTk)~ 2 at 
long limes, and the entropy and specific heat will go as 
1/((3Tk) for low enough temperatures. These results are 
in fact another manifestation of the universality prop- 
erty of this system: it implies that LL physics (with 



its ubiquitous power law dependences) cannot be mani- 
fested through the behavior of any of the thermodynam- 
ical properties, contrary to what one might expect based 
on perturbative calculations [HI ■ 

From Eqs. ([5]) and ^ we see that the dependence of 
the population becomes wider as ckfes becomes smaller 
and vice-versa, in agreement with the numerical results 
shown in Figs. [T] and [2] This has a simple interpretation: 
smaller cufes corresponds, according to the previous re- 
sults, to large g (i.e., attraction in the lead) or positive 
An. Indeed, when g is larger than 1, the local density 
of states at the edge of a LL (or at the middle of a chi- 
ral LL) diverges at the Fermi energy [31], so tunneling 
is enhanced; similarly, for An > tunneling is also en- 
hanced by Mahan's excitonic effect 17): When the level 
is empty (full) the adjacent site of the lead tends to be 
full (empty) because of the charging interaction, so tran- 
sition between these states becomes easier. In both cases, 
the population curve should indeed become broader. 

Conclusions. — To conclude, we have shown that the 
thermodynamical properties of the system are universal 
for a wide range of models, determined by only few pa- 
rameters. These properties follow a single-channel Kondo 
physics, and thus are not qualitatively affected by the 
LL phase of the lead. This implies that interesting phe- 
nomena occurring in quantum impurities coupled to LLs 
can be studied on equivalent models with non-interacting 
leads, which are much easier to study, both analytically 
and numerically (using, e.g., Wilson's numerical renor- 
malization group 25j). A clear signature of the LL phase 
can be seen when examining transport-like properties 
(e.g., the level local density of states). Alternatively, one 
could extend the model to include more than one lead. 
Both topics will be discussed elsewhere (iij ]. 
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